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A model of interacting one-dimensional fermions confined to a harmonic trap is proposed. The 
model is treated analytically to all orders of the coupling constant by a method analogous to that used 
for the Luttinger model. As a first application, the particle density is evaluated and the behavior 
of Friedel oscillations under the influence of interactions is studied. It is found that attractive 
interactions tend to suppress the Friedel oscillations while strong repulsive interactions enhance the 
Friedel oscillations significantly. The momentum distribution function and the relation of the model 
interaction to realistic pair interactions are also discussed. 



I. INTRODUCTION 



The rapid progress in cooling and magnetic trapping technology makes it conceivable to produce a neutral 

OO . ultra-cold quantum gas of quasi one-dimensional degenerate fermions following the recent achievement in three 
spatial dimensions S. Some experimental aspects of such a quasi one-dimensional fermion gas are discussed in 0. 
In three spatial dimensions, identical spin polarized fermions experience only a weak interaction because the s-wave 
' scattering contribution from short ranged pair potentials is forbidden. Scattering processes involving higher angular 
, momenta have energy thresholds and freeze out at temperatures below the Fermi temperature Tp [ jToj ] . At zero 
y—{ • temperature only a weak but long ranged dipole-dipole interaction fll| , |l^| remains. 

£N| | This classification scheme does not apply to fermions in one spatial dimension. One can only state that a contact 
interaction has no effect. It is tempting to completely neglect all interactions and treat the one-dimensional fermion 
system as an ideal gas in a harmonic trap. Some exact single particle properties of this system were given in |l3|,p|. 
The theory of Luttinger liquids p4|-^8[ (for reviews cf. ]l9|-|2l|) suggests, however, that even weak interactions destroy 
the Fermi liquid picture in one spatial dimension. This holds also true for bounded Luttinger liquids [p2[-p8[ which 
Oh, display algebraically decaying correlations. The confining potential for the cases studied so far is provided by hard 

+ii • walls. In ultra-cold quantum gases a harmonic confining potential is more appropriate. 

Besides the Luttinger approach to confined one-dimensional fermions, a wealth of exact results exist for the Calogero- 
Sutherland model [^9|,^o| with harmonic confinement and its extensions (cf., e.g., |5l]]). These results hold for the 
< 1 jr 2 pair interaction. A recent path integral approach to the thermodynamics of harmonically confined fermions ]3^ ] 
assumes a harmonic pair interaction. A different treatment of that model is given in |33| and includes results for 
Green's functions of one-dimensional fermions. The present approach, which is expected to be asymptotically correct 

K> , for large fermion numbers, is more flexible with respect to the form of the interaction. 

The method proposed is a transcription of the elementary solution of the Luttinger model (cf., e.g., [|lj) to the present 
problem. It requires, however, a different bosonization scheme which we adopt from fl34j . 

As a first application, we calculate the particle density. The particle density depends perturbationally on the in- 
teraction strength. But we are still able to present some interesting interaction effects related to the Friedel ]35| ] 
oscillations. 

The paper is organized as follows: In Sec. 2, we present the basic arguments of our approach. Sec. 3 discusses the 
occupation probabilities of single particle states and the corresponding particle density for two interaction models. 
Sec. 4 shortly discusses the momentum distribution function and its relation to the particle density. In an Appendix, 
the pertinent question of the relation between the model interaction and realistic interactions is studied and estimates 
of their respective strengths are made. 

II. LUTTINGER APPROACH 

We consider a gas of spinless identical fermions in one spatial dimension and trapped in a harmonic potential 

V(z) = \mujz 2 . (1) 
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The Hamiltonian in second quantization but without interactions is 



Hq = ^huj n c+c n , (2) 

n=0 

with one-particle energies 

Hu n = Ku)g(ri + 1/2), rt = 0,l,.... (3) 

The fermion creation operators c + and destruction operators c obey the fermionic algebra c m c^ + c^c m = S m , n . This 
ensures that each (non-degenerate) energy level e„ = huj n with (real) single particle wave function 



(4) 



is at most singly occupied. The intrinsic length scale of the system is the oscillator length I = a -1 where a is defined 
by a 2 = muji/h. H n denotes a Hermite polynomial. 

The spatial extension of the bounded system of N fermions is measured by the Fermi width L F , i.e., half the extension 
of the classically allowed region at the Fermi energy: 



L F = -y/2N-l = L„=iv-i. (5) 
a 



The Fermi energy itself is given by: 



e F = nwt(N-i) + -nu i = nwt(N-h, (6) 

with an associated Fermi wave number pc|.p| 

k F = aV2N - 1 ee ^2me F /h 2 . (7) 

The density oscillations well inside the trap can be identified |J as Friedel oscillations |H| with wave number 2k F . 
As is seen from eq. (|^), the single particle energies depend linearly on the quantum number n. A linear dispersion 
is one of the requirements which allow the treatment of the Luttinger model. Another requirement is the presence 
of an anomalous vacuum with linear dispersion extending to arbitrary negative energies. In our case, these fictitious 
free states belong to n < 0. Formally, they have wave functions tp p (z) oc T> p \^/2az) with p = —n < 0, (n: positive 
integer, T>: parabolic cylinder function) which do not belong to the Hilbert space. This deficiency is likely to have 
no effect because these functions are still normalizable in any finite spatial region [-L, L] with L ^> L F to which all 
finite range interaction effects are confined. Furthermore, one expects that for sufficiently large fermion number N 
the presence of the anomalous vacuum is harmless for processes near the Fermi energy e F oc N p8| . 
Exploiting the presence of the anomalous vacuum, it is easy to show that the density fluctuation operators 

p(m) ee Cp +m c p (8) 
p 

obey bosonic commutation relations 

\p(-p),p{q)]=p5 P ,q (9) 

for all integers p and q. The bosonic operators defined according to 



VP d+ , p > 
thus satisfy canonical commutation relations: 

[d m ,d+] = 6 m , n (11) 

for positive integers m and n. The following argument is taken over from the Luttinger model: Since the free 
Hamiltonian Hq has the same commutators with p(p): 
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[H , p{p)] = huj e pp(p), 



(12) 



as 

H = hue, md^dm (13) 

m>0 

has with the corresponding d m and d^ , the low energy (density wave) excitations of the fermionic Hamiltonian can 
as well be described in the bosonic Hilbert space. 

The form of the admissible four-fermion interaction in the full Hamiltonian is dictated by the requirement that it 
can be expressed in terms of density fluctuation operators. This restricts the form of the matrix elements of the 
interaction operator 

V = \ V(m,p;q,n) (c+c q ){c+c n ). (14) 

mnpq 

Two solvable forms were found: The matrix elements 

/oo 
dld2d3d4^ m (l)V P (2)7(l,2;3,4)^(3)^ n (4) (15) 
- OO 

must simplify according to 
a) 

V(m,p;q,n) = V a (\q - m\) <5 m _g >n _ p , (16) 
V4 = \ ^ V a {m) m {d+d m + d m d^j , (17) 



which leads to 



and 
b) 



which gives 



m>0 



V(m,p;q,n) = V b {\q - m\) 5 9 _ m ,n-p . (18) 



K = \ E ^(m) m {dl + d+ 2 } . (19) 

m>0 

V a (m) and Vf,(m) correspond to the 54 and (72 coupling functions of the Luttinger model. The Appendix discusses 

the relation of the present interactions to realistic four-fermion interactions. 

The total bosonic interaction operator, which is added to the free Hamiltonian (fL3ft, thus is 



- £ mV b (m) {d+ 2 + d r 2 „} + i mV *( m ) {dnJm + d m d+] . (20) 



V =2 

m>0 m>0 



It is diagonalized in the bosonic Hilbert space by the Bogoliubov transformation 

H = Hq + V =J2 me m f+f m + const., (21) 



m>0 

with the unitary operator 



§ = exp U S Cm - 1 • (22) 

I m>0 J 



This transformation relates the d-operators to new canonical /-operators according to 
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d m = S + fmS = fm cosh Q m - /+ sinhCm- (23) 
The transformation parameters £ m are given by 

tanh2C m = Vb }™\ , . (24) 
huje + Va(m) 

They depend also on the sign of the interactions. The spectrum of the density wave excitations is found to be 

e m = \J (huj( + V a {m)) 2 - V b 2 (m) = (%w e + V a (m j) cosh2C m - V b (m) sinh2C m . (25) 

It is also seen that V a (m) ^ only renormalizes the single particle energy hujg. For large N and at the Fermi edge, a 
local pair potential contributes equally to V a (m) and Vj,(m) because of its symmetry m <-> q or n <-> p (see Appendix). 
The matrix elements V a and Vb determine a set of coupling constants 

7 ro = sinb 2 ' 

where K m is given by 



sinh 2 Cm = ~ A [K m + — - 2 ) > 0, (26) 



hcj e + V a (m) -V b (m) 
y + Va(m) + Vb (m) 

This is in complete analogy to the Luttinger model. There is an obvious consistency condition 

\V b (m)\<\hu t +V a (m)\, (28) 



and a more subtle stability condition [fi^ 



hue + V a (m) 

We will mostly assume V a (m) — Vb(m) = V{m). Then the simple relations 

huJe v( \ 1 ( 1 



V b (m) m ^oo Q _ (2Q) 



^( m ) ^ o I — - 1 ) ( 3 °) 



hold. We will also consider two specific interaction models: 

1. A toy model, called IM1, which is defined by V(m) = V(l)(5 m ,i + <S mj _i) with a real amplitude V(l). 

2. The analogue of the usual Luttinger liquid interaction with a small exponential decay in the coupling constants 
7 m according to 7 m = exp(— r^m) 7 in order to satisfy (pf^). In the particle density considered below, further 
coupling constants a m (cf. ([40[)) appear. For these coupling constants we assume a decay according to a m = 
exp(— r a m/2) ao with ao — sgn(V(l)) 1/70 (1 + 7o)- We call this model IM2. This model is not fully consistent 
because the exponential decay of 7 m induces a non-exponential decay of a rn unless 70 -C 1 or 70 ^ 1 when 
r a — > r 7 or r a — > 2r 7 , respectively. By choosing r a properly, it is possible to restrict the error. The choice 

r 7 = r a (x V AT -1 , which corresponds to the length scale Lp, also poses no problem when N is large. 

An essential step in the applicability of the Luttinger model is the correspondence between fermionic creation and 
destruction operators and exponents of bosonic fields involving the d (or /)-operators Jl7] ^lj . The standard form of 
this bosonization is apparently not applicable to the present case. Instead we will use a method initially introduced 
in [3^ ]. The authors define an auxiliary field according to 

00 

$n{v)= e av c l =i Ja (v + 2 7 r) 7 (31) 
/— — 00 

and prove the bosonization formula for the bilinear combination 
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^(u)Mv) = G N (u - v) exp {-» (^+(u) - } exp {-* (ftu) - favj) } , (32) 

involving the non-Hermitian bosonic field 

co 1 

and the Fermi sum 

AT-1 

G w (u)= 51 e- li(u+le) . (34) 

— oo 

Equations (^), (^), j2^), and (||l]) to ([m|) provide a framework for the calculation of the fermion density and also 
of correlation functions of the density for one-dimensional fermions in a harmonic trap interacting via ( |20| ) . They do 
not, however, allow to calculate the time dependence in the single particle Green's function. 
In the remaining part we apply the concept to the fermion density in the trap. 

III. PERTURBED PARTICLE DENSITY 

The fermion density is given by the expectation value 

n(z) = $+(z)j>(z)). (35) 
Decomposing the fcrmionic destruction operator ip(z) according to 

oo 

^( z ) = ^2 ^m(z) ( 36 ) 
rn—0 

gives 

oo 

n(z) = ^2 ^ m (z)^ n (z) (c+c„). (37) 

m,n— 

Thus, the computation of the expectation values (c+c n ) is required. They become non-diagonal in the presence of 
interactions. Using the prescription of the last section, we find 

By standard methods (bosonic Wick theorem (exp [A]) = exp[(A 2 )/2] for any linear combination A of f m and /+) 
and at zero temperature (when (j m j£) = fim,n)> the expectation value ( ) = exp[— W] is calculated from 

oo 



W — W(u, v) — 2 — [y m ~ a m cosmfu + v)] {1 — costo(u — v)} , (39) 

with 



m 

m—l 



a m = sgn(V(m))^-f m (l+j m ). (40) 

In this way, a reduction to "quadratures" is achieved but nontrivial summations and integrations remain. 
W is a real and an even function of its arguments. This gives the symmetries 

(c„c„) = (c+c m ) = (c+c n )*. (41) 

It will also turn out that (c+ c„) vanishes unless \m — n\ is even. The particle density thus admits the representation: 
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oo oc M 

n(z, N,T = 0)= ^2 ^m(z) 2 (cIjCm) + 2 X! X! ^M- P [z) i>M+ P {z) (c| f _ p c M+p ). (42) 

M=0 M=lp=l 



The Z-summation in (J3^) can be performed and the two integrations appropriately transformed using the 2-k- 
periodicity of the integrand in (|38|). 

Turning to the interaction model IM1, one of the integrations can be performed, giving 
/-+ - \ 1 , 1 r j (sin((M + 1/2- N)s)\ 

exp[— 271(1 — cos(s))] 7p(2ai(l — cos(s))). 

Evidently, the factor in curly brackets, which stems from the Z-summation, provides particle-hole symmetry: For 
p =/= it has the form 

{c^N-l-M-p^N-l-U+p) = -(ctr-pCM+p), (44) 

i.e., in a plot of (c^_ p CM+p) versus M the point (N — 1/2,0) is the point of inversion symmetry. For p = one 
obtains the occupation numbers (cjjCjf), which satisfy 

{c^n-i-m^n-im) = 1 - {c^cm). (45) 

In the present fermionic case, they can be interpreted as occupation probabilities P(M) of the single particle state i\) m . 
The occupation probabilities show the expected smearing of the Fermi edge around M = N — 1 due to interactions. 
The inversion symmetry point for the probabilities is (N — 1/2, 1/2). 
Several statements concerning ( pl3| ) can be made: 

1. The result of first order many body perturbation theory in fermionic Hilbert space are exactly reproduced for 
Va^l) = by expanding ([43j) to first order in the coupling constant Vb(l) (i.e., 71 — > 0, a\ — > £1 — > Vb(l) / (2?iu>£)) 

2. The results of direct numerical diagonalization of the perturbed fermionic Hamiltonian agree with the numerical 
evaluation of ( fi"3|) up to |ao| = 1 and for particle numbers N of the order of 10, the range investigated so far 
©■ 

3. The sum rule 

dz n(z, N,T = 0)=Y / (cmCm) = N (46) 

-° M=0 

is exactly obeyed by (^): Summing the expression in curly brackets in ( [43| ) from M = to Q 3> N leads to 
the result 2nS(s)((Q +T)/2 — N). 

Returning to the particle density, the modification of P(M) and the appearance of non-diagonal contributions has 
pronounced effects: For increasingly negative values of ai the Friedel oscillations are diminished while positive ai 
increase their amplitude. This will be even more pronounced in the interaction model IM2 in which several basic 
interactions are superimposed. 

The corresponding formula for (c^ I _ p CM+ P ) is more complicated and can be given in the form: 

{o M -pCM +P ) = - S P ,o ^ [1 + Za _ cos{tro ( 47 ) 



ds f sin((M + 1/2 - N)s) \ ( Z, 



To 



2tt \ 2sin(s/2)) J \[1 + Z. ( - cos(s)] . 
[(1 + Z a - cos(i - s))(l + Z a - cos(t + s))] ao/2 . 



G 



The quantities Z 7 and Z a result from the decay functions of the coupling constants j m and a m . They are given by 

Z. ( = cosh(r 7 ) - 1, Z a = cosh(r Q /2) - 1. (48) 
Again, the sum rule (|h|) is strictly fulfilled. 

Figure 1 shows the corresponding particle density for ao = ±1, r 7 = 0.3, and r a = 0.4 obtained by numerical 



integration of ( f47| ) and by summing the series (42) up to M = 2N = 20. The occupation probabilities are smoothed 
out and more spectral weight is found above the Fermi energy. Correspondingly, the particle density outside the 
classical region is enhanced. Most startling, however, is the strong enhancement of the Friedel oscillations in case of 
strong repulsive interaction and their nearly complete suppression for strong attractive interactions. 



j n(v) 1 




FIG. 1. Dimensionless particle density n in units of the inverse of the oscillator length I versus dimensionless distance v = z/l 
from the center of the one-dimensional harmonic trap for TV = 10 interacting spinless fermions at zero temperature. Broken 
curve shows unperturbed Friedel oscillations. Strongly oscillating curve refers to a repulsive interaction with ao = 1 while the 
smooth curve is for an attractive interaction with ao = — 1. Interaction model 2 has been used. 



A discussion of the strength of the Friedel oscillations in a realistic quasi one-dimensional Fermi gas with interactions 
is addressed at the end of the Appendix. 



IV. ONE-PARTICLE MOMENTUM DISTRIBUTIONS 



Even for a confined system it makes sense to study the momentum density 

p(k) 



(CfeCfc) 



(49) 



The operator tu annihilates a fermion with (continuous) mo mentu m hk. It can be decomposed into the fcrmionic 
annihilation operators of the harmonic oscillator according to (A. 5). 

Taking into account that the expectation values (c^c n ) vanish unless \m — n\ is even, the momentum density becomes 
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P(k) = £ (/m)V5 (-l) m+ " <fe) (50) 



LXJ /ft ^ ^ 



of z — ' * — ' 

m— p— — m 



ra+p) 



Obviously, the momentum distribution always satisfies the sum rule 

dkp{k) = N. (51) 



If the expectation values (|47]) were diagonal, as in the non-interacting case, one would obtain the relation 

1 



p(k) = Arnfz = k/a 2 ). (52) 



However, the non-diagonal matrix elements are relevant: From ([43|) it is seen that in case of IM1 the matrix elements 
obey the relation (c^-_ p cjvf+ p )(— a%) — (—l) p (cj^_ p CM+p)(oix) leading to the remarkable result 

p(k, a\) = —r n(z — k/a 2 , — cti), (53) 
or 

for IM1, i.e., the momentum distribution functions behaves oppositely to the density: Attractive interactions increase 
the oscillations (with period ira 2 /kp ), repulsive interactions decrease them. This effect is also found for IM2 though 
it is less pronounced than in the density p9| ] . 

V. SUMMARY 

Using methods borrowed from the theory of the Luttinger model, we propose a scheme to treat interacting fermions 
in a one-dimensional harmonic trap analytically. As in the Luttinger model, the method rests on the linear dispersion 
of non-interacting fermion states which is extended to negative energies. A crucial step in the treatment is the 
correspondence between fermionic operators and exponentials of bosonic fields. This allows to study the fermion 
density distribution in the trap. A detailed investigation of the relation between the model interaction and realistic 
pair interactions is also given. 

Acknowledgements: The author thanks G. Alber, S. N. Artemenko, F. Gleisberg, F. Lochmann, T. Pfau, U. 
Schloder, W. Schleich, J. Voit, and C. Zimmermann for valuable discussions and Deutsche Forschungsgemeinschaft 
for financial support. 

VI. APPENDIX 

The Appendix compiles a number of formulae relating a physical two-body interaction to the four fermion matrix 
element in (|l^). We assume a local and translation invariant pair interaction 

V(zi,z 2 ; z 3 , z 4 ) = V{z\ - z 2 ) 5(zi - z 3 )S(z 2 - z 4 ), (A.l) 

and start from the representation of the interaction operator in the momentum basis: 

i r°° 

V = ^J dhdkdk'VikJc+^c+^c^Ck. (A.2) 
The Fourier transform of the interaction potential is 

/oo 
dz e lkz V(z) = V(-k) = V*(k), (A.3) 
-oo 

and the exchange symmetry of t he pa ir potential V{z) = V{—z) has been utilized. 



By rearranging the operators in (A.2), one finds 
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i r°° ~ 1 - 

V = ^ / dh dkdk' V(h) (c+ +fei c fc ) (c+_ fc » - = 0), 



(A.4) 



where iV is the fermion number operator. The self-energy contribution in (A.4), which is finite for any reasonably 
regularized potential, will be omitted. We go from the momentum basis to the harmonic oscillator basis by means of 



Cfe 



n=0 



with H 



1 



dze ikz ip n (z) = —ip n (z = -j). 



/27T J-oo ' a ' a 1 

Thus the usual form for the matrix elements, as used in the interaction operator 



1 °° 

^=9 ^Kk?. 11 ) fe)(c p + c„), 



mnpq—O 



namely 



V{m,p\ g, n) 



dz\dz2 ipm(zi)ip p {z2)V(zi - Z 2 )lpq(zi)lp n (z2) 



(with the obvious symmetries m <-> q, n «-> p) takes on the less familiar form 



V(m,p;q,n) = (-!)»+« -L 



dfc , V'(fci){/*- kl /*}{/*' +fcl /5'} 



V(m,p;g,n) = (-l) n_N i m+n +f+« 



2tt 



*i C'„ l< j(A: 1 )y(fc 1 )Cp„(-fci), 



with (q > m is assumed) 



a* 



k-h\ ( k 



(A.5) 



(A.6) 



(A.7) 



(A. 



(A.9) 



It can be further evaluated by performing the convolution in (A.9). By using the method described in fl9j, one obtains 

1 



(A.10) 



(All) 



I ml ( k\ 



q—m 



exp 



Aa 2 



ql \aV2 

In order to eliminate a possible contact interaction in V(z), we set 

V cS (k) = V(k) - V(k = oo) 



4r m) (^) = (-ir +q c qm {k 1 ). 



(A.12) 



to obtain the following exact expression for the matrix elements of the potential in the basis of oscillator states and 
for q > m, n > p: 



V{m,p; q, n) = i^+v+i + (-1)""™}, 



' m\p\ a 
2n\q\ 2tt 



(A.13) 



p OO 

/ dv e ~ V«- m +«-P- 1 )/ 2 V eS (k = ay/2v) L^ m ^ (v) (v). 
Jo 

These matrix elements are of particular importance for values of the indices m, n, etc., which correspond to wave 
functions near the Fermi level, i.e., m ~ N etc.. Setting q — m + Q and exploiting m 3> Q (e.g., Q < V~N), the 
asymptotic expansion Q] of the Laguerre polynomials can be applied: 
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L<£Hv)~^Y Jq(2V^)- (A.14) 
This formula implies the approximation (m + Q)l ~ rrftm\. Setting n = p + R, one finds 



V(m, p;m + Q,p + R)~ iQ+*U3-+±Vl |L 



(A.15) 



jf e-* V eS (k = aV2^j Jq(2V^) J R {2y/pv). 



The integration region near v — will contribute little to matrix elements with Q > and i? > (the Q = R = 
case in ( |l6| , |i"8"| ) is irrelevant) provided V(k) is bounded for small wave numbers k. In order to estimate these matrix 
elements, we can thus assume v to be not less than of order unity and apply the asymptotic expansion of the Bessel 
functions J. This gives: 

V(m,p;m + Q,p+R)^iQ +R t^+S-^l^ [°° ^ e ~4 (A.16) 



V2 7T 2 J a U 



cos ( uy/m - ^ 7T - '- ) V eS ( k = ^ ) COS ( M^P 



where a is of order unity. We choose a = 1. Finally, setting 7? = Q 4- AJ?, leads to 

V(m,p;m + Q,p + Q + AR)^ — - cos(— -) / - e -TM(ifc = -= ) (A.17) 



tt 2 V2 2 J x u V2 

{cos (u(Vto - y/p)) + sin (u(s/m + y/p)) } 



Equation (A.17) admits the following discussion of the matrix elements: 

1. The matrix elements depend weakly on m and p in the ranges |Am| = \m — N\ < \/N and |Ap| = \p— N\ < \^N. 



2. Elements with AR = are (nearly) independent of Q < v N and all have the same sign: They are positive for 
repulsive interactions and negative in the opposite case. 

3. Varying AR, produces alternating signs which weakens the effect of these matrix elements. 

These are essentially the properties which justify the model interactions (|l6| ) and (|l8|). In this connection, it should be 
mentioned that, as in Luttinger liquid theory, one does not need a faithful representation of the physical interaction 
|jl8| , flb"f to capture essential features of the interaction. 

For a numerical evaluation of the interaction strength, we begin with the interaction of longitudinally aligned magnetic 
dipolcs (cf. jll],[l2]]) in a quasi one-dimensional magnetic trap for which the effective (attractive) potential is 

"<*> = -f" 2 o?+W <A ' 18) 

Here, an ultra-violet cut off has been introduced. It is given by the width of the one-dimensional channel, i.e., 
d = a^ 1 , where at is the reciprocal harmonic oscillator length of the transverse (ground state) wave function. This 
length is large on an atomic scale. Introducing the ratio A = u>i/u>t, one has d = y/X/a. In case of a completely filled 
one-dimensional trap, which we he ncefor th assume, N = 1/A holds. 
Fourier transforming the potential ( A. 18 ), gives 

V(k) = V (V\-)K 1 (VX-). (A.19) 
a a 

This expression involves a modified Bessel function and is bounded (though not analytic) at k — 0. It vanishes for 
large momenta. The prefactor contains the magnetic moment fi of the fcrmions and reads: 
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V 



Mo P 2 ^ 2 



For IM1, the quantity 
1/(1) ; 



/°° X X u 

<My - Kx{uyj -)} exp{- — } cos (u(y/m- v /p)) 



(A.20) 



(A.21) 



is needed. The contribution from the second term in ( A. 17 ) containing yfm + y/p is negligible. The factor in square 
brackets can be written as 



V 



1 

71 



hue 



(A.22) 



where ?7M is the atom mass. The prefactor of the scaling energy fwji is estimated for 6 Li to be —3 • 10 3 (A^/10 4 ) 
adopting the data from p). The remaining integral in (A.21) is about unity for |m — p\ < \/~N. So it needs about 
N = 10 6 atoms in the trap to produce a significant effect of the dipole-dipole interaction. 

For scattering processes in three dimensions, the van der Waals potential is more relevant. The dilute atoms in the 
trap (mean distance larger than 100 Bohr radii) explore mostly the long-range part of it. Reducing the R~ 6 potential 
to an effective potential in one dimension, gives: 



with the Fourier transform 



V(z) = 



7T A 



(z 2 + d 2 ) 3 ' 



Following p 



V(k) = -^e~ kd (3 + 3kd + k 2 d 2 ). 
the coefficient A can be expressed in terms of the coefficient Cq therein, giving 



V(l) 



c 6 



16tt\/2 



2 5 5 

e a a t a 



huj£ 



du 



exp(-it 2 /4- uy/X/2) 



(A.23) 



(A.24) 



(A.25) 



(3 + 3uv/a72 + m 2 A/2) cos(u(v^- y/p)) , 

where eg and a are electron charge and Bohr radius, respectively. Fixing ujg at 2n ■ 10s _1 , as in |j, the evaluation of 
( A.25| ) gives a total prefactor of hu>i according to w —6 • 10 -7 (A^/10 4 ) 5 / 2 . Even for N = 10 6 , this is small compared 
to the dipole-dipole interaction. This reversal in the importance of the two interactions can be traced back to the 
ultra-violet cutoff d = a^ 1 needed in the quasi one-dimensional forms of the interactions. 

The choice of transverse cut-off determines the strength of the effective one-dimensional interaction. Our choice 
is motivated by the following fact: The transverse trap direction is associated with only one intrinsic length scale 
(Lp)t = l/a t — h — IvX ~ l/kp and this length coincides with the shortest intrinsic length scale of the longitudinal 
trap direction. 

The amplitude of the Fricdcl oscillations inside a bounded Fermi sea scales with the particle number N as 1/N. 
Adopting the numbers estimated above, it is impossible to detect Friedel oscillations in a gas of 10 6 6 Li atoms with 
present techniques. However, the interaction strength V(l) between magnetic dipoles scales as f^m 3 ^ 2 . Using 53 CV 
instead of 6 Li enhances the interaction by nearly three orders of magnitude. Fermionic molecules are even more 
promising, especially when they are polar (cf., e.g., (4^]): An electric dipole moment of 1 Debye (p = aoeo/2.5) leads 
to an interaction strength which is enhanced by a factor (0.8/as) 2 « 10 5 (qj: fine structure constant) over that from 
magnetic dipoles of magnitude /is {^b- Bohr's magneton). 

Thus it seems that the issue of Friedel oscillations in an interacting quasi one-dimensional degenerate Fermi gas is 
not purely academic. 
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